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ABSTRACT 

We present the thermal evolution and emergent spectra of solidifying terrestrial planets along with 
the formation of steam atmospheres. The lifetime of a magma ocean and its spectra through a steam 
atmosphere depends on the orbital distance of the planet from the host star. For a type-I planet, 
which is formed beyond a certain critical distance from the host star, the thermal emission declines 
on a timescale shorter than approximately 10® years. Therefore, young stars should be targets when 
searching for molten planets in this orbital region. In contrast, a type-H planet, which is formed inside 
the critical distance, will emit significant thermal radiation from near-infrared atmospheric windows 
during the entire lifetime of the magma ocean. The Ks and L bands will be favorable for future direct 
imaging because the planet-to-star contrasts of these bands are higher than approximately 10“^-10“®. 
Our model predicts that, in the type-H orbital region, molten planets would be present over the main 
sequence of the G-type host star if the initial bulk content of water exceeds approximately 1 wt%. In 
visible atmospheric windows, the contrasts of the thermal emission drop below 10“^® in less than 10® 
years, whereas those of the reflected light remain 10“^® for both types of planets. Since the contrast 
level is comparable to those of reflected light from Earth-sized planets in the habitable zone, the 
visible reflected light from molten planets also provides a promising target for direct imaging with 
future ground- and space-based telescopes. 

Subject headings: infrared: planetary systems, planets and satellites: terrestrial planets, radiative 
transfer 


1. INTRODUCTION 

Theoretical studies predict that the final stage of ter¬ 
restrial planet formation involves a series of giant impacts 
between protoplanets. This stage is triggered by the dis¬ 
persal of nebula ga s, and probably lasts approximatel y 
10 to 100 Myr fe.g. lAgnor et al.lll999t IChambersll2001ll . 
In the aftermath of the last im pact, the pla nets proba¬ 
bly begin in a molten state ('e.g. lCanugl2004ll . Detecting 
molten terrestrial planets in extrasolar systems is of great 
significance in testing the widely accepted view of their 
hot origins. 

The rocky planets can acquire water throughout for¬ 
mation by various processes, such a s the local accretion 
of hydrated and icy p l anetesimals (iMuralidharan et al.l 
|20 o 4 iKing et ^120101 : iMachida fc Abell2010D . c hemical 
reactions between nebula g as and silicate melts (ISasakil 
llQQOtITkoma fc Gendall200CT . and delive ry via ice-bearing 
bodies from the outer orb i tal region (IMorbidelli et al.l 
120001: [Raymond et al.l[2007t iWalsh et al.ll201lh . A suR 
stantial amount of water would be delivered during plan¬ 
etary formation. Water acquired through the lattermost 
process alone would exceed the current Earth’s ocean 
mass by more than tens times. Since gi ant impacts do 
not blow off the entire water inventory (jGenda fc Abel 
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1200511 , early molten planets are likely to be covered by a 
steam-rich atmosphere. 

Emergent spectra of hot molten planets have been in¬ 
vestigated for several atmosph eric compositions, includ¬ 
ing water-rich atmospheres. iMiller-Ricci et abl (|2009f) 
first presented thermal emission spectra for the hot 
surface conditions, considering various possible atmo¬ 
spheric compositions; solar and enhanced metallicity at¬ 
mospheres, H2O-CO2 atmospheres, and Venus’ compo¬ 
sition atmospheres. Since the solidus temperature of 
peridotite rocks is approximately 1,400 K at the plan¬ 
etary surface (e.g. iTakahashU 1 1 98611 . the molten surface 
emits strong radiation not only in the infrared (IR), but 
also in the visible and near-IR. They showed that, al¬ 
though the steam-rich atmosphere obscures the hot sur¬ 
face, intense thermal emission from the surface would 
leak through near-IR at mospheric window s, which are 
relatively transparent. iLupu et al.l (|20I4ll further ex¬ 
plored emergent high-resolution spectra for atmospheric 
compositions in equilibrium with magmas with compo¬ 
sitions of Earth’s continental crust and the bulk silicate 
Earth, taking into account atmospheric chemistry and 
structure self-consistently. They found out that the at¬ 
mospheres are dominated by H2O and CO2, considering 
quenching effects. Additional opacity sources, such as 
CH4 and NH3, could reduce the brightness temperature 








































2 


by a factor of two to three, but the emergent spectra still 
reveal the features si milar to those with out these addi¬ 
tional opacity sources. iLupu et all (|2014ll concluded that 
the most favorable wavelength range for direct imaging 
with future telescopes will be 1 to 4 /rm. 

These previous studies revealed that the mass of the 
steam atmosphere is a determining factor for spectral de¬ 
tectability of hot molten planets. However, the mass of 
the steam atmosphere does not remain constant during 
the solidification process. The atmosphere would grow 
by degassing fr om the planetary interior or would es¬ 
cape into space (lElkins-Tantonll201lt ILebrun et al.ll2013l : 
iHamano et al.l 120131) . Atmospheric evolution would sig¬ 

nificantly affect the emergent spectra during solidifica¬ 
tion of the magma ocean. 

In addition, atmospheric evolution is also crucial to 
the th er mal evolution of magma oceans (iM atsujfcAbd 
1986allbl: iZahiile et a 1 119881 lE lkins-Tantonl 120081 l2niic 
Lebrun et al.l l20l'^ Hamano et al.l 1201311 . Since the 
strong greenhouse effect of water vapor prevents heat 
from escaping, a magma ocean is sustained longer un¬ 
der a more massive atmosphere. Although a growing 
atmosphere keeps the planet in a molten state for a 
longer period of time, its spectral signatures become 
fainter as a result of increasing surface obscuration. Con¬ 
versely, if the atmosphere escapes with time, the planet 
becomes brighter, but the magma ocean quickly solidi¬ 
fies over a shorter timescale. Thus, a trade-off relation¬ 
ship is expected between the spectral intensity and oc- 
cu rrence rate of the rnagma- covered planets, as suggested 
bv iMiller-Ricci et'all (1200911 . 

Variations in water inventory would also affect the 
planetary spectra and their time evolution. The water 
in the current Earth’s ocean amounts to approximately 
0.023 wt % in bulk content, and its total, including the 
water stored in t he crust and mantle, is probably a t most 
0.23 wt % (e.g. lOhtanil [20051 : iHirschiriannl l2006ll . Ex¬ 
trasolar terrestrial planets could exhibit large variations 
in water amount just after formation as a result of a 
combination of water supply and loss processes. Recent 
calculations have indicated that the lifetime of magma 
oceans depends on the initial wat er inventory and th e 
orbital distance from the host star (|Hamano et al.l[2013ll . 
Therefore, the occurrence rate of hot molten planets as 
a function of orbital distance would be related to the ex¬ 
pected variety in initial water endowment on terrestrial 
planets. 

As described above, spectral calculations of hot molten 
planets have been presented for various atmospheric com¬ 
positions by sev eral authors (|Miller-Ricci et akl 120091 : 
ILudu et ^1201411 . They have concluded that the de¬ 
tectability in terms of spectral intensities depends pri¬ 
marily on the mass of the atmosphere. On the other 
hand, recent evolutionary calculations have shown that 
the mass of the steam atmosphere overlying a magma 
ocean gre atly varies with time by degassing and escape 
processes (|Hamano et al.l[2M^ . They also indicated that 
the evolutionary history of the steam atmosphere and 
its cooling time strongly depends on the orbital distance 
from the star. Although they suggested that the occur¬ 
rence rate of molten planets depends on the planetary 
orbit as well, their evolutionary model is not able to pre¬ 
dict the variations of thermal emission spectra along with 
the evolution of the steam atmosphere because of the as¬ 


sumption of a gray gas in their radiative transfer model. 

In the present paper, we incorporate a newly developed 
non-gray radiative trans f er mo del into the evolutionary 
model bv IHamano et al.l (j2013ll to describe the spectral 
evolution of hot molten planets along with the evolu¬ 
tion of steam atmospheres during the solidification of a 
magma ocean. This enables us to self-consistently exam¬ 
ine the brightness variations of solidifying planets and 
the duration of the molten state. We also investigate 
the dependence of the brightness of solidifying planets 
on orbital distance and obtain a consistent relationship 
between planet-to-star contrast and orbital separation. 
In this paper, we consider Earth-sized terrestrial planets 
around a Sun-like star. 

In Section[2l we summarize our proposed model and its 
parameters. We describe the spectral features of thermal 
emission and the albedo of molten planets with a steam 
atmosphere. In Section|3l we describe the time evolution 
of surface temperature, atmospheric pressure, and ther¬ 
mal emission spectra. In Section |4l we present diagrams 
of planet-to-star contrast and orbital distance for near- 
IR atmospheric windows and discuss future prospects for 
direct imaging. We also address the possibility of putting 
constraints on water inventory originating from the for¬ 
mation stages of terrestrial planets from the predicted 
occurrence rate. Color variations of the solidifying plan¬ 
ets are presented as well. In Section [5l we compare our 
model spectra with those obtained by previous models 
and discuss parameter uncertainties. In Section 0 we 
summarize our main findings and present our conclu¬ 
sions. 

2 . MODEL DESCRIPTION AND SPECTRAL 
FEATURES 

2.1. Radiative transfer calculations 

The evolutionary model used in the p resent study 
is that described in IHamano et al.l (|2013f) . except for 
the calculations of radiative transfer. In the present 
study, we replace a g ray radiative transfer code used in 
IHamano et al.l (l2013fl with a newly developed, non-gray 
code to calculate thermal emission and reflection spec¬ 
tra from the planet. In the non-gray code, a two-stream 
formulation is used to calculate the radiative transfer in 
a plane-parallel atmosphere. We calculate planetary and 
stellar radiation separately in a wavenumber range from 
0 to 30,000 cm“^ with a spectral resolution of O.OI cm“^. 
The albedo of the planetary surface is fixed at 0.2. See 
the Appendix for details on the radiative transfer model 
and its validity. 

A steam atmosphere is divided into 100 layers, which 
are equally spaced on a logarithmic scale of pressure. 
The pressure at the top of the atmosphere is fixed at 
I Pa in all calculations. The temperature structure is 
prescribed such that it follows an adiabat ic lapse rate in 
all atm ospheric layers, based on results of lAbe fc Matsuil 
(|1988f) . Eor atmospheric layers in which the temperature 
is above the critical temperature of water or the rela¬ 
tive humidity is less than I, the lapse rate is assumed 
to be a dry adiabat. For the other layers, the lapse rate 
is given by a pseudo-moist adiabat under the assump¬ 
tion that condensates are removed immediately as they 
form. For simplicity, the effects of clouds are ignored in 
the present study. Since we consider a pure steam atmo- 
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Fig. 1.— (a) Thermal emission spectrum through a steam atmosphere of 50 bar with a surface temperature of 2,500 K. The light-pink line 
indicates a high-resolution spectrum with a wavenumber resolution of 0.01 cm“^, corresponding to a resolving power Rx of approximately 
1,000,000 at 1 fxm. The red line indicates the average thermal radiation such that Rx becomes 100. Gray bands represent atmospheric 
windows, which are defined in Table ^ (b) High-resolution (pink) and average (red) spectral albedo for a stellar zenith angle 6^ of 
cos”^(l/\/3). (c) Pressure level at which the optical depth is unity for absorption (red line) and for scattering (dotted gray line), (d) 
Atmospheric structure. 


sphere, the pseudo-moist adiabat overlaps the saturation 
curve of water in the present study. We apply a liquid- 
vapor saturation curve at temperatures above 273 K and 
an ice-vapor curve at temperatures below 273 K. The dry 
adiabat is calculated using th e Peng-Robinson equ ation 
of state in a manner similar to I Abe fc Matsuil (1198811 . but 
the heat ca pacity of water vapor in the ideal gas state is 
taken from lWagner fc Prufl (|2002ll as a function of tem¬ 
perature in order to make the heat capacity applicable 

to higher-temperature co nditions. __ 

We use HITEMP 2010 (iRothman et alllM^ as spec¬ 
tral database to calculate the absorption coefficients 
for lines of water vapor. For all calculations, a Voigt 
line profile is used to consider the combined effect of 
Doppler and pressure b r oaden ing, using a computa¬ 
tional code by iHumlfcekI (jl982l l. Although no reliable 
model for water continuum is available for the high- 
temperature, high-pressure c onditions considered herein, 
we adopt a MT_CKD 2.5 (iMlawer et al.l 1201211 model 
for its continuum together with a wavenumber trunca¬ 
tion of 25 cm“^, which is similar to a procedure used in 
recent studies on radiative tra nsfer calculations of hot 
and thick steam atm ospheres (iKopparapu et al.l 120131 : 
iGoldblatt et al.l[2013ll . As the number of line transitions 
in the HITEMP 2010 exceeds 10®, we generate a grid of 
absorption coefficients of water vapor and calculate the 
optical depth of each layer by interpolating the absorp¬ 
tion coefficients. 


We consider Rayleigh scattering by water vapor with 
a scattering cross-section per molecule of 


327r® 1 / mg — 1 

3 V ^0 


2 


n 5 


( 1 ) 


assuming that the polarizability volume of water vapor 
does not depend on temperature or pressure. In the 
above equation, Wn is the wavenumber, Na is the Avo- 
gadro constant, and rig and mg are the molar volume 
and the refractive index in the reference state, respec- 
tively . We adopt a refractive index of 1.000254 from 
lAllenI (II973I1 for D lines of Na at the standard tempera¬ 
ture and pressure (0°C, 1 atm). 


2.2. Thermal emission 

Figure[IJa) shows the thermal emission spectrum from 
0.33 to 100 ^m with a surface pressure of 50 bar and a 
surface temperature of 2,500 K. At this surface temper¬ 
ature, the peak wavelengths of the Planck function B\ 
and \B\ are approximately 1.2 and 2 /im, respectively. 
The planetary surface emits strong radiation in the vis¬ 
ible and near-IR. Although the hot planetary surface 
emits significant radiation, the outgoing thermal emis¬ 
sion from the top of the atmosphere is greatly attenu¬ 
ated by absorption and scattering by the steam atmo¬ 
sphere. The stronger thermal radiation passes through 
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TABLE 1 

The band parameters used in the present study. 


Band 

Ac [^m] 

AA [^m] 

u 

0.365 

0.068 

B 

0.44 

0.098 

V 

0.55 

0.089 

R 

0.70 

0.22 

I 

0.85 

0.18 

Y 

1.04 

0.20 

J 

1.26 

0.24 

H 

1.65 

0.29 

Ks 

2.16 

0.32 

L 

3.55 

0.57 

M 

4.77 

0.45 

N 

10.47 

5.19 


References. — The band parameters for I, Y, and 
J are defined by au thors, and the other band par ame- 
ters are taken from ITraub &: Qppenheimerl Il2010l h 


the wavelength regions corresponding to the atmospheric 
windows, at which the absorption of water vapor is rela¬ 
tively weak. 

Figure djc) shows the atmospheric pressures at which 
the optical depth is unity for absorption and scattering, 
respectively. In the far-IR and absorption bands of water 
vapor, the upper layers, in which the pressure is typically 
less than 10 “^ bar, are responsible for the thermal emis¬ 
sion spectrum, whereas the atmospheric windows reveal 
the deeper and hotter layers. At wavelengths longward 
of approximately 0.6 fim, the pressure level at which the 
optical depth of absorption is unity is much lower than 
that of scattering. In the near-IR wavelength region, 
therefore, the brightness temperature provides a good 
estimate of the atmospheric temperature at an optical 
depth of unity. 

At the visible wavelength, the scattering process has a 
dominant role in attenuating the thermal emission. At 
shorter wavelengths, the deepest layer could contribute 
to the emergent spectrum because of the absence of 
strong absorption lines. The actual probed depths, how¬ 
ever, depend on the atmospheric amount and continuum 
absorption of water vapor, the physics of which remain 
poorly understood. If the atmosphere is thin enough to 
allow the thermal radiation emitted from the surface to 
pass through the atmosphere, a Rayleigh scattering slope 
appears in the planetary spectrum at sufficiently opti¬ 
cally thin wavelengths. Measuring the scattering slope 
can, in principle, provide information about the atmo¬ 
spheric column mass, the surface temperature and the 
surface albedo. 


2.3. Spectral albedo 

Figure [TJb) shows the spectral albedo for a stellar 
zenith angle 6 ** of cos“^ (l/v^- The spectral albedo 
is defined as the ratio of the upward diffusive radiation 
to the downward direct radiation at the top of the at¬ 
mosphere. The bond albedo of the planet is obtained as 
a weighted mean of the spectral albedo by the Planck 
function with an effective stellar temperature. 

In addition to the thermal radiation, the spectral 
albedo is higher within the atmospheric windows. More¬ 
over, the spectral albedo varies greatly at the visible 
wavelengths, from approximately 1 at 0.4 /rm to less 
than 0.1 beyond 0.9 pm. Therefore, in the near-IR, the 


reflected stellar light is expected to have a smaller contri¬ 
bution to the emergent spectrum. Although the spectral 
and resulting Bond albedos increase with stellar zenith 
angle, the difference is less than 0.2 throughout the vis¬ 
ible wavelengths between 0* of 0 and 80 degrees. In 
the present study, we use the spectral and Bond albe¬ 
dos with 0 * of cos“^ (I/a/S) degrees as a representative 
value. Spectra composed of both reflections and emis¬ 
sions are presented later herein under conditions that are 
consistent with the atmospheric evolution in the course 
of solidification of the magma ocean. 


2.4. Evolutionary model 

Energy and water budgets during plan etary solidifica¬ 
tion a re treated in the same manner as iHamano et al.l 
dM^. Here, we briefly summarize our proposed model, 
focusing especially on the major change s made in the 
presen t study. For additional details, see IHamano et al.l 
dMl). We assume an adiabatic temperature profile 
in the planetary interior to obtain heat capacity of the 
whole silicate portion. The heat flux from the magma 
ocean is calculated as the difference between the outgoing 
planetary radiation and the net incoming stellar radia¬ 
tion. In order to reduce the computing time for the evolu¬ 
tion calculations, the planetary radiation and albedo are 
precalculated on 882 surface-temperature-pressure grid 
points, from 1,350 to 3,000 K and from 5 x 10“^ to 
5,000 bar, respectively. In order to follow the evolution of 
planet-to-star contrast, the thermal radiation from each 
band is also tabulated after averaging the high-resolution 
spectra with a resolving power of 100 , which is a typical 
value for future direct imaging. We consider a Sun-like 
host star, the bolometric luminosity and XUV radiation 
of which evolve in a manner simila r to solar anal ogs. We 
adopt a solar standard model by iGoughI (|198ir ) for lu¬ 
minosity enhancement with stellar age and linearly ex¬ 
trapolate it after a stellar age of 4.567 Gyr. The stellar 
spectrum is assumed to be a blackbody spectrum with 
an effective stellar temperature of 5,800 K. 

Water is assumed to be partitioned between the at¬ 
mosphere and the magma ocean according to its solubil¬ 
ity into basaltic melts. The incorporation of water into 
cumulates is also considered according to the partition 
coefficients of water between melts and solids and by as¬ 
suming 1% of trapped interstitial melts. The mass loss 
rate of water is calcu lated from an energy -limited escape 
flux of hydrogen Fh (|Watson et al.lll9^ with a heating 
efficiency 77 of 0 . 1 , as follows: 


Fh = 47ri?p[ 


Rp\ 

GMpi 




-XUV 


4 



( 2 ) 


where i?pi and Mpi are the planetary radius and mass, 
respectively, a is the orbital distance from the host star 
in AU, and G the gravitational constant. Here, 
is the XUV flux from the host star at 1 AU. The stel¬ 
lar XUV flux decreases rapidly with time, whereas the 
ratio of the X-ray luminosity to the bolometric lumi¬ 
nosity remai ns nearly constant at 10 ~^ for young 
active stars (IVilhu fc Waited Il987t iStauffer et al.l 11994 
iPizzolato et al.ll2003f) . Based on a power-law r elatio n de- 
rived from o bservational data by iRibas etldl (1200511 and 
[RibasI (l2009ll . we adopt the following expression for the 
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Fig. 2.— Evolution of a planet at (a) 1 AU (type I) and (b) 0.7 AU (type II) with an initial water mass of 5 M^q. The third and bottom 
panels, respectively, show the thermal radiation and the maximum value of the planet-to-star contrast within each atmospheric window. 
The dotted vertical line in each figure indicates the time at which the magma ocean has solidified, which is approximately 3.9 Myr for the 
type-I planet and 100 Myr for the type-II planet. 


evolution of 5”^^^ as a function of stellar age r in Gyr: 

/ / \ - 1.23 


^xuv 


^xuv / ^ \ ^ 

/ \ - 1.23 

j for r > 


( 3 ) 


and 


r,at = l- 66 x 102 OL-O- 64 , (4) 

where Lboi is the bolometric luminosity of the host star in 
erg s^r^ and Tq and 5'™^, which are the age and XUV 
flux at 1 AU of the current Sun, are fixed to 4.567 Gyr 
and 4.59x10“^ J 8“^ m“^, respectively. The value of Tsat 
is approximately 0.066 Gyr for the solar standard model 
by IGoughI (jlQSlf) . The mass loss rate of water is nine 
times as large as the mass loss rate of hydrogen under 
the assumption that all the dissociated oxygen atoms 
are consumed to oxidize the abundant su rface magma 
(|Gillmann et al.l[2009l : iHamano et al.ll^OlSf) . 

We consider terrestrial planets with the same mass and 
bulk composition as Earth. We start the evolution calcu¬ 
lations from a very hot state with a surface temperature 
of 3,000 K and stop the calculations when the surface 
temperature reaches the solidus temperature at the sur¬ 
face (approximately 1,370 K) or the stellar age reaches 
10 Gyr, which is comparable to the main sequence life¬ 
time of the host star. The solidification is assumed to 
start at a stellar age of 0.05 Gyr (tq). We consider the 
initial total inventory of water ranging from 0.01 to 50 
times the current ocean mass of Earth Meo of 1.4 x 10^^ 
kg. This corresponds to a range in bulk water content of 
from 2.3 x 10“"^ to 1.2 wt %. Furthermore, we consider 


the orbital distance from the star of between 0.4 and 1.5 
AU. 


3. SPECTRUM EVOLUTION DURING 
SOLIDIFICATION 

A hot planet cools and solidifies as it emits large out¬ 
going radiation. Spectral variations during planetary so¬ 
lidification reflect the evolution of surface temperature 
and atmospheric pressure. At early stages of evolution, 
the highest fraction of water is expected to dissolve in 
a deep magma ocean b ecause of the high solubility of 
water into silicate melts (|Zahnle et al.ll201^ . The steam 
atmosphere therefore accounts for a modest fraction of 
its total inventory. As the solidification proceeds, the sil¬ 
icate magma becomes enriched in water and degases at a 
rate that d epends on the solidifi cation rate of the magma 
ocean fe.g. lElkins-Tantonll2008ll . Whether the steam at¬ 
mosphere grows with time depends on the net balance 
of the degassing rate of water fro m the interior and it s 
loss rate by hydrodynamic escape (jHamano et al.l[^013H . 
Conside ring a magma ocean c overed with a steam atmo¬ 
sphere, IHamano et al.l (j2013l l proposed that terrestrial 
planets have a distinct early evolution that depends on 
whether the orbital distance from the host star exceeds a 
critical distance. The critical distance OciEI is defined as 
the distance at which the net stellar heating flux is equal 
to the radiation limit Fn^ of steam atmospheres: 


Gcr 


0.83 


( film 

V280Wm-2y 1^0.7 5'^°' 


1/2 


U-0.2y 


45) 


^ We rescaled the expression for acr in IHamano et al.l (120131) 
using the typical values in the present study. 
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(a) 1 AU (b) 0.7 AU 




Fig. 3.— Spectral evolution in the visible and near-IR during the course of solidification of the planet at (a) 1 AU (type I) and (b) at 0.7 
AU (type II) with an initial water mass of 5 Mgo • The color of the thermal emission spectra represents the elapsed time after solidification 
starts. The orange line indicates the thermal spectrum at the time at which the planet has solidified. 


where Ofpi is the Bond albedo of the planet, and 5^°^ and 

are the bolometric fluxes at 1 AU of the host star 
and the current Sun, respectively. 

When the troposphere is saturated by water va¬ 
por to a sufficiently deep level, the outgoing plan¬ 
etary flux becomes constant. This constant value 
is referred to as the radiation limit and corresponds 
to the lower limit of the planetary radiation for the 
hot surface conditions considered in the present study. 
The radiation limit of the steam atmospheres has 
been extensively investigated, a nd its value i s ap¬ 
proximately 2 8 0 to 310 W m~ ^ (I Abe fc Matsuil 119881: 
Kastind 119881: iNakaiima et al.l 119921 : iKopparapu et al l 
2013t iGoldblatt et al.l 120131) . Our model y i elds a radi¬ 
ation limit of 280 Wm“^. iHamano et al.l (|2013r) also 
classified terrestrial planets into two types according to 
their orbital distance: type I for planets beyond Ucr and 
type II for planets inside Ocr- These two types of planets 
will solidify over vastly different timescales, and their at¬ 
mospheres also evolve differently, due to different cooling 
mechanisms. In the following subsections, we describe 
the spectral changes over the solidification of type-I and 
type-II planets with an initial water inventory of 5 Meo 
as a nominal case. 


3.1. Type I: planets beyond Ocr 

The type-I planet is characterized by a sufficiently far 
distance from the host star, where the net incident stel¬ 
lar radiation is smaller than the radiation limit. Since 
the outgoing planetary radiation never falls below the 
radiation limit, the type-I planet is self-luminous over its 
solidification period. The minimum heat flux, which is 
given by the difference between the radiation limit and 
the net incident stellar flux, determines the overall solid¬ 
ification timescale, which is approximately 4 Myr in this 
nominal case (Fig. HJa)). The rapid solidification re¬ 
sults in the monotonic growth of the steam atmosphere 
towards the end of the magma ocean period by keeping 
the degassing rate of water higher than its loss rate. 

Although the surface temperature during the magma 
ocean period is extremely high, the bolometric luminos¬ 
ity is not so large compared to that from the hot surface. 
Here, we define an effective emission temperature Tgff, 


using the planetary radiation Fpi, as follows: 

Teff = ^Fpi/a. (6) 

The effective emission temperature never exceeds 800 K 
in this nominal case and is typically as low as 300 K 
over the solidification period. The thick steam atmo¬ 
sphere thus sustains the magma ocean longer, whereas 
the planet itself looks as faint as a planet with a cool or 
moderate surface temperature. 

The most recognizable diagnostic sign of a hot surface 
is emission at the visible and near-IR wavelengths. Al¬ 
though the massive steam atmosphere strongly absorbs 
the thermal emission from hotter and deeper layers, ther¬ 
mal emission can leak through atmospheric windows. 
The third panel of Fig. [2](a) shows the time variations 
of thermal radiation from atmospheric windows. In the 
early stages (less than 10^ years), the thermal radiation 
from all bands exceeds 10^ Wm“^/rm“^. These high 
fluxes, however, do not last long. At most of the at¬ 
mospheric windows, the thermal fluxes rapidly decrease 
with time and fall below 10“^ Wm“^/xm“^ due to the 
decrease in the surface temperature and the growth of 
the steam atmosphere. The thermal radiation of the M, 
L, and Ks bands initially decreases with time as well, but 
eventually levels off. 

Figure [31(a) shows the spectral evolution in the visible 
and near-IR. The planetary thermal radiation declines 
with time, and then, at wavelengths longer than approx¬ 
imately 2 /im, the thermal radiation becomes constant. 
This is a result of the same mechanism as the occur¬ 
rence of the radiation limit of the steam atmospheres 
and indicates that, for the wavelengths at which water 
vapor has strong absorption, the emission level becomes 
as high as that at which the thermal structure is uniquely 
determined by the saturation curve of water. Since the 
thermal structure around the emission level remains the 
same, the thermal radiation from the M, L, and Ks bands 
becomes constant, as shown in Fig. |2](a). 

3.2. Type II: planets inside Ocr 

In contrast to the type-I planet, the type-II planet re¬ 
ceives a net stellar radiation that is larger than the radi¬ 
ation limit because of its proximity to the host star. The 
planet therefore would reach a radiative energy balance 
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during the course of solidification. Early in its solidifi¬ 
cation, the outgoing planetary flux far exceeds the net 
stellar flux, and the planet is highly self-luminous, as in 
the case of the type-I planet (Fig. EKb)). In the case of 
the type-II planet, however, the planetary flux decreases 
with time and eventually balances the net stellar flux 
at approximately 0.6 Myr. If no water loss occurs, the 
radiative balance is perfectly achieved, and the magma 
ocean does not solidify. Actually, as the steam atmo¬ 
sphere becomes thinner with time due to the hydrody¬ 
namic escape, the planet emits slightly larger radiation 
than the net stellar radiation. The type-II planet thus 
slowly cools and solidifies, losing its atmosphere. 

The solidification rate of the magma ocean is regulated 
by the loss rate of water such that a sufficient net loss 
of the steam atmosphere occurs so that the net outgoing 
flux becomes positive. The resulting solidification time 
is approximately 100 Myr in this nominal case, which is 
much longer than that of the type-I planet. The most 
striking feature is that the higher band fluxes continue 
over this long solidification period (third panel in Fig. 
[2Kb)). The thermal radiation from the atmospheric win¬ 
dows ceases to decrease at ~ 0.6 Myr, at which radiative 
energy balance has been roughly achieved. The subse¬ 
quent spectra do not vary with time at the opaque IR 
wavelengths or the atmospheric windows, which are rel¬ 
atively transparent (Fig. EKb)). 

The long-lived and constant thermal emission in the 
window regions occurs by attaining quasi-energy balance 
on the type-II planet. In the quasi-energy balance state, 
the outgoing planetary radiation roughly balances the 
incoming net stellar radiation. In order to satisfy this 
requirement, the amount of the atmosphere must be reg¬ 
ulated during the course of solidification. As the surface 
temperature decreases with time, the amount of the at¬ 
mosphere must decrease such that the planet maintains 
a planetary radiation that is nearly equal to, or, strictly 
speaking, slightly larger than, the net stellar flux (Fig. 
HKb)). Consequently, the resulting spectra do not change 
during a quasi-energy balance state, except at the shorter 
wavelengths, at which the opacity of the atmosphere is 
thin enough for the outgoing thermal emission to be sen¬ 
sitive to surface temperature (Fig. |3Kb)). 

The above requirement for planetary radiation gener¬ 
ally applies to planets in (quasi-)energy balance, irre¬ 
spective of the atmospheric amounts and composition. 
Although additional opacity sources would increase the 
opacities at some wavelengths and would change the 
overall spectra, the planet in quasi-energy balance is still 
expected to emit high thermal radiation from the rela¬ 
tively transparent wavelengths. 

Figure |3Ka) shows the evolutionary paths of surface 
temperature and pressure (the Ps — path) during so¬ 
lidification for various initial inventories of water. The 
atmospheric pressure increases at early stages due to the 
high degassing rate. For an initial water inventory larger 
than 0.2 to 0.5 Meo, tbe quasi-energy balance is attained 
at some particular surface temperature. The Ps—Ts paths 
in the quasi-energy balance states are approximately in¬ 
dependent of the initial water endowment, as long as the 
planets solidify during the main sequence of the host star. 

For a larger initial inventory of water, the planet 
reaches quasi-energy balance at a higher surface temper¬ 
ature and solidification requires a longer period of time. 



Surface Temperature [K] 
Initial Bulk Content of Water [wt%] 



Fig. 4.— (a) ps — Tg evolutionary paths of the planet at 0.7 AU. 
The number on each curve indicates the initial inventory of water 
in the current Earth’s ocean mass Meo of 1.4 x 10^^ kg. The thick 
parts indicate that the planet is in quasi-energy balance. Here, the 
‘quasi-energy balance state’ is defined as the state in which the 
relative difference between the planetary and net stellar radiation 
is less than 5%. The hooks of the ps — Tg curves with an initial 
water inventory of 20 and 50 Meo result from the enhancement 
of the stellar luminosity, (b) Solidification time of the planet at 
0.7 AU as a function of initial water inventory. The filled circles 
represent the cases in which the quasi-energy balance state has 
been attained at some point during solidification, and the open 
circles indicate the cases in which the quasi-energy balance state 
has not been attained at some point during solidification. The gray 
dashed line is the time required for complete loss of water. 


In particular, for the case in which the quasi-energy bal¬ 
ance has been achieved during the course of solidification, 
the lifetime of the magma ocean is well approximated by 
the time required for total loss of the initial water inven¬ 
tory (Fig. EKb)). If the planet is endowed with water 
equivalent to more than 16 MeO; the molten state would 
exist for over 10 Gyr, which is comparable to the lifetime 
of the G-type host star. In this case, the surface temper¬ 
ature gradually increases with time in accordance with 
the luminosity enhancement of the star. 


4. DETECTABILITY OF MAGMA-COVERED 
PLANETS WITH A STEAM ATMOSPHERE 


Dedicated space coronagraphic missions, such a s 


the Terrestrial Planet Finder C_(|Levine_et_^ 


Occulting Ozone Observatory ( Savranskv et al 


2009), 


2010 1. 


WFIRST-AFTA (ISoergel et all 120151) 7 and Exo-C and 
Exo-®, have been conducted over long periods of time. 


^ Exo-C and Exo-S Final Reports in 2015, 
http: //exep.jpl.nasa.gov/stdt/ 














Required Minimum Telescope Diameter [m] (d=10 [pc]) Required Minimum Telescope Diameter [m] (d=10 [pc]) 



Required Minimum Telescope Diameter [m] (d=10 [pc]) Required Minimum Telescope Diameter [m] (d=10 [pc]) 




Elapsted Time I 

10'* 10^ 10® 10^ 10® [year] 

Fig. 5. — Contrast ranges of the thermal emission during the magma-ocean period at each orbital distance with an initial inventory of 
water of 5 Mbq: (a-) J, (b) H, (c) Ks, and (d) L bands. Each color bar indicates a contrast range during solidification, with the tick 
marks indicating the elapsed time for each order of magnitude. The cross symbols indicate the contrasts at which the magma ocean has 
completely solidified. The numerical value under each color bar in Fig. ISfa) is the overall solidification time. At orbital distance smaller 
than 0.8 AU, the planet-to-star contrast increases with time at the very last stage of solidification because of the low degassing rate. 
The increase in the contrast at the very last stage is neglected in this figure due to its short duration relative to the total solidification 
time. The upper horizontal axis shows the telescope diameter i^tel 3-t which the angular separation is equal to the diffraction limit, 
0dl ~ 0.01 X (Ac/0.5 ^m) X (iltei/10 m)“l^, for the case in which the distance to the system d is 10 pc. 


The primary goal of such missions is to detect re¬ 
flected light in visible bands from Earth-sized or super- 
Earth-sized planets around solar-type stars. These mis¬ 
sions will be capable of detecting the reflected light 
from an Earth-sized planet that is typically 10^° times 
fainter than the light of the host star. Moreover, 
recent progress in extreme adaptive optics has im¬ 
proved the limiting contrast of direct imaging using 
CTOund-based telescope s, such as Ge mini Planet Image r 
([Macintosh et al.l 120141) , SPHERE (iZurlo et al.l p0141 . 
Palomar Project 1640 dOppenheimer et al.l 120121) . and 
SCExAO (|Jovanovic et al.l 1201 3D . These ground-based 
high-contrast instruments have an advantage of the small 
inner working angle on the future 30- to 40-m class tele¬ 
scope. These missions will open the possibility of detect¬ 
ing molten exoplanets. In this section, we discuss the fea¬ 
sibility of the direct detection of molten planets based on 
the capabilities of these future missions. We also discuss 
characterization via direct imaging, which would provide 
evidence of the hot origins of molten planets as well as 
the basis for more extensive spectroscopic studies for the 
purpose of characterization of these planets. 


We consider two important quantities that character¬ 
ize the ease of detection of a planet against stellar speck¬ 
les. One is planet-to-star contrast, which is defined by 
the ratio of the planet flux to the stellar flux, and the 
other is the planet-star angular separation on the celes¬ 
tial sphere. The latter can be converted to the semi¬ 
major axis a for a given distance assuming an observer- 
star-planet angle (phase angle) of 90°. For future space 
coronagraphic missions, we regard a contrast of 10“^*^ at 
0.1 arcseconds of angular separation, which corresponds 
to an Earth-twin at 10 pc, as a benchmark for detectabil¬ 
ity. For future ground-based high-contrast instruments, 
although it is difficult to set a benchmark, we fiducially 
regard 10“^ — 10“® of the contrast at 0.01 to 0.03 arc- 
seconds (corresponding to the diffraction limit of the J 
- L bands for a 30-m telescope), as a benchmark be¬ 
cause one of the main goals of these missions is to de¬ 
tect Earth- or super-Earth-sized planet s in the habitable 
zone around late-type stars (se e e.g. iMacintosh et ahl 
2006t iMatsuo fc Tamural 120101: iKawahara et al.l 20121 : 


Guvon et al.ll2012HCrossfieldll2013HMales et al.ll2014lb 
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Required Minimum Telescope Diameter [m] (d=10 [pc]) Required Minimum Telescope Diameter [m] (d=10 [pc]) 



Orbital Distance [AU] 


Orbital Distance [AU] 


Required Minimum Telescope Diameter [m] (d=10 [pc]) 


Required Minimum Telescope Diameter [m] (d=10 [pc]) 



Orbital Distance [AU] 


Orbital Distance [AU] 


Planet-Star Contrast _ i 

^ 0 ''^ 10'" 10’’° tO'^ 10'® 10'^ 10'® 

Fig. 6. — Contrast evolution of the thermal emission at each orbital distance with an initial inventory of water of 5 M^o- (n) Jj (b) H, 
(c) Ks, and (d) L bands. Color shows a planet-star contrast during solidification. The upper horizontal axis in each panel is the same as 
Figure [5] 


4.1. Flux contrast and orbital separation 

Direct detection requires that the planets have a suffi¬ 
ciently high planet-to-star contrast. As described in the 
previous section, the time evolution of the thermal emis¬ 
sion spectra is closely related to the atmospheric evolu¬ 
tion during solidification, which strongly depends on the 
orbital distance from the star. Figure [5] shows the con¬ 
trast evolution of the thermal radiation in near-IR win¬ 
dows for various orbital distances in the case of an initial 
water inventory of 5 Meo- Although the total contrast 
consists of the thermal radiation and the reflected stellar 
light, in most cases, the thermal radiation is dominant in 
the near-IR band, especially for type-II planets 1 114.3|) . 

The contrast evolution of type-I planets is different 
from that of type-II planets. For type-I planets, which 
are located beyond Ocr, the range of the planet-to-star 
contrast is approximately independent of the semi-major 
axis. The planet-to-star contrast is initially as high as 
I0“® to I0“®, and then decreases on a common timescale 
of less than ~ 10® years (Fig. [0]). This is because the 
magma-ocean lifetime of type-I planets is determined by 
the minimum heat flux, which is approximately indepen¬ 
dent of the orbital distance. In the case of a smaller 
inventory of water, the planet-to-star contrast would in¬ 
crease as the atmosphere becomes thinner, whereas the 


timescale for decreasing planet-to-star contrast becomes 
shorter. In the case of a larger inventory of water, 
the thermal radiation reaches the radiation limit earlier. 
This prolongs the duration of the magma ocean by a fac¬ 
tor of at most 2. 

The contrasts of type-II planets first decrease with time 
as well, but then cease to decrease during the course 
of solidification (Fig. [SJ. The minimum value of the 
contrast of type-II planets is higher than that of type-I 
planets. Since the planets are in quasi-radiative energy 
balance, they spend most of the solidification period at 
the minimum contrast level in this nominal case (Fig. |6|). 
The minimum contrast level is insensitive to the initial 
water inventory. As long as the planet solidifies on a 
shorter timescale than the main sequence lifetime of the 
host star, the minimum contrast level has the highest 
likelihood of detection. For the case in which the magma 
ocean is sustained for as long as the main sequence of 
the star, the planet-to-star contrast gradually increases 
with time from the minimum contrast level, because the 
planet must emit a larger planetary radiation from the 
near-IR band in order to maintain quasi-energy balance 
in accordance with the luminosity enhancement. 

The minimum contrast increases for the planet clos¬ 
est to the star, because this planet receives a larger net 
stellar flux and so must emit larger planetary radiation. 
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Consequently, planets closer to a star are more favorable 
for direct imaging in terms of spectral contrast. On the 
other hand, the planet closest to the star also receives 
a stronger stellar XUV flux, which fuels hydrodynamic 
escape. Since the solidification time of type-II planets is 
approximately equal to the loss time of water, the occur¬ 
rence rate of molten terrestrial planets is expected to be 
higher at orbital distances closer to Ocr- 

As shown in Figures and m the type-II molten 
planets maintain relatively high contrasts in the near-IR 
band. The type-II planets with a < 0.5 AU are favorable 
targets for future ground-based direct imaging, because 
the contrasts of the Ks and L bands are above the levels 
of ^ I0“® and ^ 10”^, respectively, for the entire life 
of the molten state. The duration with the contrast > 
10 “® is approximated with the overall solidification time, 
which becomes longer with the initial water inventory, 
as shown in Section 14.21 Although these planets have 
smaller angular distances, in the near future, a 30- to 40- 
m ground-based telescope will resolve them within tens 
of parsecs. On the other hand, the near-IR contrast of 
the type-I planets declines below ^ 10“® and ^ 10“^ 
on a short timescale of ^ 10^ years for the Ks and L 
bands, respectively. Young stars should be targeted when 
searching for type-I planets. 

The thermal radiation in the near-IR band overwhelms 
the reflected stellar radiation for type-II planets. The 
contrast of the emission in the visible bands rapidly de¬ 
creases to the 10“^° level by a few 10^ and ^ 10® years for 
the type-I and type-II planets, respectively (see Figure 
[2D. Therefore, with respect to the visible band, the re¬ 
flected light should be considered in order to discuss the 
total cont rast of the planets. These issues are discussed 
in Section 1231 

4.2. Relation between predicted occurrence rate and 
water endowment 

For the type-II planets, the lifetime of the magma 
ocean strongly depends on how much water is acquired 
during planet formation. This suggests the possibility of 
placing constraints on the initial water inventory based 
on the observed occurrence rate of hot molten planets. 
Figure [7] shows an orbital region where molten Earth¬ 
sized planets potentially exist as a function of stellar 
age. The predicted lifetime of the magma-covered plan¬ 
ets shows a sharp transition around Ocr- The type-I plan¬ 
ets during the magma ocean period have less probability 
of being detected, whereas the probability of detecting 
molten type-II planets is greatly enhanced with initial 
water inventory. 

The calculated lifetimes of type-I and type-II planets 
allow us to estimate expected number of G-type stars 
which host Earth-sized molten planets. For G-type stars 
whose distance from the solar system is less than the 
distance d, the expected number of G-type stars with 
molten planets fVoCd) can be estimated as follows, 

A+tot 

iVD(d)^iVG(d)-PEG-^^, (7) 

where NG{d) and fc^d) are, respectively, the number 
and average age of G-type stars within the distance d 
from the solar system. Peg is the fraction of G-type stars 
with terrestrial planets in a range of planetary orbital 


Initial bulk content of water 
1 _ 2 5 _ 10 _ 20 [Meq] 


0.023 0.047 0.12 0.23 0.47 [wt%] 



0.4 0.6 0.8 1 1.2 1.4 1.6 

Orbital Distance [AU] 



0.4 0.6 0.8 1 1.2 1.4 1.6 

Orbital Distance [AU] 


Fig. 7.— Orbital region where molten planets potentially exist 
as a function of stellar age in a range from 0.01 to 10 Gyr (a), with 
the assumption that the solidification starts at stellar age tq of 
0.05 Gyr. The panel (b) is a zoom-in of the panel (a) for a stellar 
age younger than 0.1 Gyr. Each color represents the initial bulk 
content of water required for a molten state to be sustained for a 
given stellar age and orbital distance. As well as tq, the age of the 
Sun Tq, 4.567 Gyr, is indicated by an arrow on the right side. The 
critical distance with planetary albedo api of 0.2 is also shown in 
each panel. 

distances considered, is the total duration during 

which terrestrial planets are in a molten state. Let us 
consider the case when d is 10 pc. The total number of 
G-type stars is then around 30. The average age of G- 
type stars is a s sumed to be 5 Gyr. Using Kepler data, 
IPetigura et al.l (I2013f) reported that the fraction of GK- 
type stars which host Earth-sized planets is ~ 0.12 for 
orbital distances of 0.05-0.42 AU, and ^ 0.086 in the 
habitable zone. Here, we adopt Peg of 0.1 for both type-I 
and -II planets. 

Type-I planets solidify on a timescale of the order of 
10 ® years, which shorter than typical time intervals of 
giant impacts at late stage of formation. In this case, 
the total duration would be expressed as a product of an 
average number of giant impacts nci and the duration 
of a molten state after one collision Atmoi- According to 
N-body simulations, Earth-like planet s would experience 
10-20 giant impacts during forma tion (jKokubo & Gend^ 
I2010I l^ewart & LeinhardtJ 12012( 1. We adopt riQi of 10. 
Using a typical duration of a magma ocean of 2 x 10® 
years, we obtain ~ 0.01 for Nu{d = 10pc) in the type-I 
orbital region. This value is common for type-I planets 
which formed with initial water inventory comparable to 
or exceeding the total amount of water of the modern 
Earth. The potentially favorable target for molten type- 
I planets would be young and other spectral type stars 
such as Fomalhaut. On the other hand, type II plan¬ 
ets have a wide range of solidification time after single 
impact, which strongly depends on the initial amount of 
water. In the extreme case that the solidification time 
exceeds the period of giant impact phase, which probably 
lasts until the stellar age of 0.1 Gyr, the total duration of 
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a molten state is given by the duration of a molten state 
after single collision. In the type-II orbital region, the 
initial amount of water thus greatly affects the probabil¬ 
ity of detecting molten planets. Given that does 

not exceed 10 Gyr, Nuld = 10 pc) ranges from ~ 0.01 
to 6 in the case of type-II planets, considering the initial 
inventory of water of from one-tenth to several tens times 
the Earth’s ocean mass. 

Here, we address the uncertainties in parameters which 
affect the loss time of water, the timing of the last giant 
impact, the heating efficiency and the EUV flux from the 
G-type host star. The giant impact phase is triggered 
by the dispersal of nebula gas and lasts until a stellar 
age of approximately 0.1 Gyr. Therefore, the last gi¬ 
ant impact would occur between the stellar ages of 0.01 
and 0.1 Gyr. The heating efficiency 77 in Eq. ([2]) de¬ 
pends on the structure and details of the energy budget 
in the upper atmosphere and remains poorly constrained. 
The estimated values ha ve a large spread and typically 


range from 0.1 to 0.5 le.g.lWatson et al.lll981l: 

Chassefierd 

ll99filKulikov et al.ll2007l: IShematovich et al.l 

201411. The 


saturation level of XUV flux and its ending time Tgat 
remain quite uncertain. If we simply assume that the 
XUV-to-bolometric-luminosity ratio scales as the X-ray- 
to-bolometric-luminosity ratio, which is observationally 
estimated to be io-3-2± 0-3 during the satura tion phase 
for one solar-mass star (|Pizzolato et al.l [200311 . the vari¬ 
ation of the XUV saturation level is nothing more than 
a factor of 3 and Tsat ranges from 0.06 to 0.13 Gyr. Al¬ 
though the XUV flux could exhibit greater fluctuation 
for pre-main sequence stars, its effects might not be so 
significant because of its short duration relative to the 
subsequent evolution. Although no observational data 
are available at wavelengths between 360 and 921 A due 
to strong absorption by the interstellar medium, the to¬ 
tal XUV flux would not have profound uncertainty com¬ 
pared with that of the other parameters, because the es¬ 
timated en ergy flux in this w avelength range is at most 
10 to 20% (jRibas et al.ll2005( l. 

Although our model cannot predict the exact inventory 
of primordial water due to the numerous uncertainties de¬ 
scribed above, estimating the order of magnitude of the 
exact inventory of primordial water would be still use¬ 
ful. Our model suggests that, if the initial bulk content 
of water exceeds approximately 1 wt%, molten planets 
would be present over the main sequence of the host G- 
type star (Fig. [7]). Since the critical distance, which 
separates the type-I and type-II planets, originates from 
the presence of the radiation limit of steam atmospheres, 
detection of a critical distance would offer a diagnostic 
test to determine whether water is generally one of the 
major volatile species of extrasolar terrestrial planets. In 
the orbital region of type-II planets, the detection of hot 
signatures in the planetary spectrum will provide a lower 
limit for initial water inventory, which is sufficient to sus¬ 
tain the molten state for a given stellar age, whereas no 
detection will constrain its upper limit, so that the planet 
will lose most of its water and become solidified. 

4.3. Reflected light and color variations of solidifying 

planets 

The actual spectrum from the planet is obtained as the 
superposition of its thermal emission and reflected stellar 
light. The thermal radiation generally overwhelms the 
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Fig. 8.— Superimposed spectra and band color of a planet at 
0.6 AU with an atmospheric pressure of 100 bar and a surface 
temperature of 2,000 K. Under this condition, the planet is in quasi¬ 
energy balance. The thin black lines in the top panel indicate the 
planet-to-star contrast of superimposed spectra. The spectrum 
components are indicated by the red line for the thermal emission 
and the green line for the reflected light, respectively, for a phase 
angle of 90°. The middle panel shows the color index of each band 
to the V band relative to the host star. As a reference, color indexes 
relative to the Sun ar e plotted in the bottom panel for all planets 
(ILundock et al.l 120091) . The superscript attached to the planet’s 
name denotes the observational period of the data used: 1 for May 
2008 and 2 for 2008 Nov. 

reflected stellar radiation during the highly self-luminous 
stage. Figure in shows the planet-to-star contrast of the 
superimposed spectra of the planet in quasi-energy bal¬ 
ance. According to the Lambert assumption, the planet- 
star contrast Cpf{X) of the reflected light is given as 
follows: 

C”'(A) = ^a(A)(:^)" (8) 

= (9) 

TT 

where a is the spectral Bond albedo, and is the 
Lambert phase function with respect to the phase angle 
/3 = Z(star — planet — observer). The overall superposed 
spectrum in the quasi-energy balance state is V-shaped 
with a bottom around I /rm, except for the phase angle 
very close to the new phase (Fig. |8]). The dominant com¬ 
ponent is thermal emission at wavelengths longer than ^ 
I fxm and is reflected stellar light at the shorter wave¬ 
lengths. This originates from the difference in the peak 
wavelength of the Planck function between the surface 
temperature and the stellar effective temperature. The 
















































12 



H-K Color 

Fig. 9.— Color-magnitude relation of the thermal emission from 
the solidifying planets with a steam atmosphere at 1 AU (green) 
and 0.7 AU (red) with an initial water inventory of 5 Meo- The 
dashed lines indicate data for the cases in which the stellar light 
reflected by clouds is taken into account with a phase angle of 90°. 
Here, we assume that the clouds are uniformly distributed and the 
coverage and albedo are 0.5. We adopt the H-K color of the Sun 
(0.05) as that of the host star. The brown circles are observed 
brown dwarfs IjLeggett et al.ll2002l : IKnaop et ajjl2 004lb The letters 
represent the planets i n the solar system l|Lundock et al.l 1200911 . 
The results reported bv ILupu et al.l 1120141) are also shown for total 
surface pressures of 10 bar (filled triangles) and 100 bar (open 
triangles) for surface temperatures of 2,200 and 1,600 K for the 
case in which the atmospheric composition is in equilibrium with 
magma having the composition of the Earth’s continental crust. 

peak wavelengths are approximately 1.4 fim for 2,000 K 
and 500 nm for 5,800 K. 

For P = 90°, the contrast of the reflection can be esti¬ 
mated as follows: 

where i?© is the Earth radius. Because the albedo in the 
U-B band is ^ 1, as shown in Figure [Tib) and the reflec¬ 
tion dominates in this band, we roughly expect contrasts 
of ^ 10“^° and ^ 10“® for the type-I and type-II plan¬ 
ets, respectively. These contrasts are within the scope of 
the future space direct imaging survey. The contrasts de¬ 
cline steeply from the visible wavelengths to the near-IR 
wavelengths, which reflects the deep atmosphere of the 
planets. 

The reflected light is bluer than that of the host star, 
whereas the thermal emission in the near-IR is redder 
than that of the host star. The middle panel in Fig. [8] 
shows the color index of each band to the V band rela¬ 
tive to the host star as a function of phase angle. As the 
phase angle increases, the color becomes redder for all 
bands because the contribution of the thermal emission 
increases. The color indexes of the I and R bands re¬ 
main bluer than that of the host star, whereas the color 
indexes of the other spectral bands becomes redder than 
that of the host star. The color indexes could vary by 
approximately 6 for the J, H, Ks, and L bands during a 
complete cycle of orbital motion. These color changes are 


comparable to the color differences between rocky plan¬ 
ets, such as Earth and Mars, and gaseous planets such 
as Uranus a nd Neptune. 

Following ILudu et all (1201411 . we consider the color- 
magnitude diagram of the molten planets. Figure [5] 
shows the evolutionary paths of the thermal emission in 
the H-K color and absolute H magnitude diagram. The 
molten planets are as faint as the planets in the solar 
system, but are distinctly redder in color. This is be¬ 
cause the thermal emission, which is much redder than 
the stellar light, dominates the planetary spectrum for 
molten planets, whereas the color of the planets in the 
solar system mostly reflects the reflected light. We also 
estimate the effects of stellar light reflected by clouds 
on the H-K color by assuming that the cloud coverage is 
0.5. In this case, as the thermal emission becomes fainter, 
the reflected light becomes more dominant in the near- 
IR spectrum. The reflected light could determine the 
color at the later stage for the planet at 1 AU, whereas 
the color is still much redder than that of the star dur¬ 
ing solidification of the planet at 0.7 AU, including the 
quasi-energy balance state. A self-consistent evolution¬ 
ary calculation, including cloud formation, will be the 
subject of future research to examine actual color varia¬ 
tions that occur during the course of solidification. 

5. DISCUSSIONS 

We compare the bright ness temperatures o f our m odel 
with those report ed by iMiller-Ricci et al.l (12009(1 and 
ILudu et all (|2014ll for wavelengths ranging from 1 to 
5 ^m (Fig. Iini). Although the atmospheric composi¬ 
tions differ among thes e models, our mod el spectrum is 
more similar to that in ILudu et al.l (|2014ll . especially for 
the atmospheric composition in a chemical equilibrium 
with silicate melts having the composition of the Earth’s 
continental crust. The brightness temperature, i.e., the 
thermal radia tion, of the our m odel tends to be lower 
than that of iLupu et al.l (120141) for continental-crust- 
composition magmas, probably because of the H 20 -rich 
composition and the lower temperature in the upper at¬ 
mosphere of our model. 

Although we have thus far discussed only pure steam 
atmospheres, the inclusion of other gaseous species would 
not substantially alter our conclusions. Adding other 
opacity sources would reduce the planet-to-star contrasts 
in atmospheric windows for type-I and -H planets in 
a self-luminous state, whereas doing so would only af¬ 
fect the relative contrast of atmospheric windows for the 
type-H planets in quasi-energy balance. If the thermal 
radiation decreases at certain wavelengths during the 
quasi-energy balance state, the planet must emit larger 
thermal radiation at other wavelengths in order to main¬ 
tain the energy balance while reducing the amount of 
the steam atmosphere. The same would be said of the 
effects of clouds. At very early stage, a molten planet 
would have no cloud cover because of the hot temper¬ 
ature. As the temperature drops and the atmosphere 
grows, however, clouds would start to form by conden¬ 
sation of water vapor in the atmosphere. Clouds would 
increase the planetary albedo, which makes the critical 
distance smaller by reducing the net incident stellar ra¬ 
diation. Clouds also would enhance a greenhouse effect, 
which makes the critical distance larger. The presence of 
clouds therefore would affect the actual position of the 
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Fig. 10.— Brightness temperatures obtained by the different 
models. The total atmospheric pressure is 10 bar in all models, 
where as the surface te mperature i s 1,500 K for IMiller- Ricci et alJ 
1120091) and 1,600 K for ILupu et al.l II2014I) and our model. We con¬ 
sider a pure steam atmosph ere. The model atmosphe r e consists of 
90% H 2 O and 10% CO 2 for IMiller- Ricci et al.l 1120091 ). ILupu et al.l 
lf20TD considered atmospheres in equilibrium with the hot surface 
having the composition of the Earth’s continental crust and bulk 
silicate Earth (BSE). Although these model atmospheres are com¬ 
posed by various gase o us sp ecies, they consist primarily of H 2 O 
and CO 2 . ILupu et al.l II2Q14) stated that, by taking into account 
quenching effects, the spectra for BSE composition magmas be¬ 
come similar to that of the Earth’s continental crust. 

critical distance through a balance of the opposing effects 
of cloud albedo and a greenhouse effect. However, as long 
as a planet is type II, it would emit strong thermal radi¬ 
ation through the atmospheric windows even with some 
cloud cover, to maintain energy balance for the same rea¬ 
son as the case with additional gaseous opacity described 
above. 

To a first order approximation, we have assumed that 
the planetary interior has an adiabatic temperature pro¬ 
file in our evolutionary model. Although horizontal in¬ 
homogeneity would exist in the convective magma ocean, 
especially around boundary layers, the planet would 
have a nearly uniform surface temperature under a thick 
and massive steam atmosphere considered in the present 
study, because the timescale of heat transfer by atmo¬ 
spheric circulation is generally much shorter than the 
timescale of radiative cooling and of heat transfer by con¬ 
vection of the magma ocean. 

In the evolutionary model, we have also assumed that 
the magma ocean effectively transports internal heat 
to the surface and that the atmosphere is in solution 
equilibrium with the magma ocean for calculations of 
a degassing rate. Since the assumption of the solution 
equilibrium leads to a maximum estimate of the atmo¬ 
spheric growth rate, the thermal radiation predicted in 
the present study would be a minimum estimate, whereas 
the decline timescale of the emergent spectra would be 
overestimated. 

Heat flux and a degassing rate from a magma ocean 
strongly depends on the magma viscosity. Thermal 
structure and heat tran sfer of a magma o cean has been 
extensively reviewed by iSolomat^ (j2007tl . Viscosity of 
ultramafic silicate melt near liqnidu s is as low as 0.1 Pa 
s at low pressures (e.g. iShawlIlQTlIl . The similar value 
can be assumed for the magma viscosity under the man¬ 
tle conditions because of small temperature and pres¬ 
sure effects for a low-viscosity and completely d epolymer- 
ized melt (see discussions bv iSolomatovI (|20n7[ ll. As the 
Rayleigh numbers are typically in the range of 10^°-10^° 


with the low viscosity (|Solomatovl l2007t lElkins-TantonI 
1200811 , the magma ocean at early stage is expected to be 
highly turbulent and in a state of vigorous convection. 
The highly turbulent magma ocean is considered to be 
homogeneous and its thermal structure is to be nearly 
adiabatic. 

Along with cooling, the magma ocean starts to solid¬ 
ify at great depth, because the slopes of liquidus and 
solidus of Earth’s mantle are steeper than thos e of adia- 
bats in the magma ocean (lAndranlt et al.ir20II|l . On the 
other hand, when a solid crust starts to form depends on 
the energy balance at the surface. The energy balance 
between heat flux from the magma ocean and radiative 
heat flux from the surface determines the surface tem¬ 
perature. A steam atmosphere has a strong blanketing 
and greenhouse effect. Even with a surface atmospheric 
pressure of 10 bar, its total opacity exceeds 10® in the 
10-/rm atmospheric window. Therefore, the insulating 
atmosphere maintains the surface in a molten state if 
the heat flux from the magma ocean is sufficiently high. 

The most important parameter for heat transport of 
the magma ocean is the magma viscosity, which varies 
dramatically depending on the melt fraction in the 
magma. Experimental studies have shown that viscosity 
of partial melts abruptly increases from a liquid-like low 
value to a solid-like high value (rheological transition), 
when its melt fractio n falls below a critical value, approx - 
imately 40 % (e.g. Ivan der Molen and Patersoa 1197911 . 
As long as the surface temperature is high enough to 
keep the melt fraction larger than the critical value, the 
magma ocean can transport the internal heat to the sur¬ 
face efficiently by vigorous convection, and also keep a 
high renewal rate of the surface. As the surface temper¬ 
ature decreases, the melt fraction decreases and reaches 
the critical value, and then solid-state convection takes 
over. The heat flux from the magma ocean is no longer 
high enough to keep the surface molten. The surface tem¬ 
perature rapidly drops a nd a conductive l i d crust starts 
to de velop at the surface (|Solomatovll2007t iLebrun et al.l 

[2?fn . 

As the solidification proceeds, the depth of the magma 
ocean becomes shallow and the total mass of the magma 
decreases. The magma ocean becomes enriched in water, 
which is incompatible to silicate cumulates. The solu¬ 
tion equilibrium of water between the atmosphere and 
the magma ocean is assumed in our model, whereas ac¬ 
tually a certain degree of supersaturation of water might 
be require d for degassing by bub bling and H2O diffusion 
processes. iHamano et al.l (|2013[1 discussed the effects by 
the simple treatment for the degassing processes on the 
overall solidification time of the magma ocean. For de¬ 
gassing involving the bubble formation, bubbles are re¬ 
quired to grow fast and detach from the magma flows. 
Although the supersaturation required for the efficient 
bubble detachment remains uncertain, even using 50-fold 
water solubility shortens the solidification time by a fac¬ 
tor of at most 4, based on the results by IHamano et al.l 
(j20I3[l . They also estimated the supersaturation required 
for the replenishment of the atmosphere against its loss 
by the molecular diffusion process only, usi ng a molecu¬ 
lar dif fusion coefficient in basaltic melt by iZhang et al.l 
(|2007ll . They showed that a modest degree of supersat¬ 
uration is enough for the replenishment, as long as the 
magma viscosity is a liquid-like low value. Eddy diffusion 
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process would significantly lower the required supersat¬ 
uration. However, at a very late stage after the melt 
fraction reaches the critical value, both the resurfacing 
rate and the efficiency of bubble detachment from the 
magma flows would be greatly reduced due to the solid¬ 
like high viscosity. Redefining the solidification time as 
the time required for the surface temperature to drop 
to 1,700 K, around which the melt fraction is estimated 
to reach the critical fraction for the rheological transi¬ 
tion, shortens the lifetime of the magma ocean on type-1 
planets by a factor of less than 5 in the parameter range 
considered herein. This estimate is consistent with the 
typical du ration of the magma ocean phase of 1 Myr re¬ 
ported by iLebrun et al.l (|2013ll . For type-11 planets, the 
redefinition of the solidification time has less effect on 
the lifetime of the magma ocean, as long as the water 
endowment is sufficiently large so that the quasi-energy 
balance is attained at a surface temperature above 1,700 
K. These uncertainties would be smaller compared to 
the parameter uncertainties associated with the energy- 
limited escape rate, as discussed above. 

Although we have examined only Earth-sized planets 
in the present study, the orbital-dependent dichotomy 
in lifetime and spectral evoluti on would exist for differ - 
ent planetary masses as well. iKopparapu et al.l (|2014[1 
reported that the radiation limit increases by approxi¬ 
mately 20 % for a rocky planet with 5 times the Earth’s 
mass. The critical distance therefore would become 
closer to the host star for larger planetary masses, com¬ 
pared to that for Earth-sized planets. Although the 
larger radiation limit also would have a role to make 
the cooling of a type-I planet faster, the overall dura¬ 
tion of a molten state may be roughly proportional to 
the planetary mass because of the larger heat capacity 
of the planet (jElkins-Tantonl[2011tl . On the other hand, 
for a massive type II planet, the lifetime of a magma 
ocean would become extremely longer, with the assump¬ 
tion that the bulk content of water is the same. This is 
because the total amount of water is proportional to the 
planetary mass, whereas the total escape rate of hydro¬ 
gen is independent on the planetary mass itself, as long as 
the energy-limited escape rate is considered. In terms of 
spectral detectability, more massive planets are expected 
to be brighter and there fore easier to detect beca use of 
their larger surface area (|Miller- Ricci et al] 1200911 . Fu¬ 
ture studies will expand this line of research by taking 
into account planetary mass, different types of stars and 
an expected variety of atmospheric compositions. 

6 . CONCLUSIONS 

We have examined the lifetime of the magma ocean and 
its spectral evolution for Earth-sized terrestrial planets 
along with the evolution of the steam atmosphere. Due 
to the presence of the radiation limit of steam atmo¬ 
spheres, the thermal and spectral evolution of the magma 
ocean strongly depends on the planetary orbital distance 
from the host star. Planets beyond the critical distance 
flcr (type I) are self-luminous throughout solidification. 
The thermal radiation from atmospheric windows de¬ 
creases on a timescale shorter than ^ 10® years. The 
lifetime of the magma ocean is less sensitive to the semi¬ 
major axis in the type-I orbital region. Consequently, 
young stars should be targets in searching for molten 
planets in this orbital region. 


In the case of a planet formed inside Ocr (type II), the 
outgoing planetary radiation can balance the net incom¬ 
ing stellar radiation after a short self-luminous phase. 
The emergent spectra during the quasi-energy balance 
state are approximately independent of the surface tem¬ 
perature and pressure. The type-II planet maintains 
high thermal radiation from near-IR atmospheric win¬ 
dows throughout the lifetime of the magma ocean. The 
lifetime of the magma ocean increases for a larger initial 
inventory of water. 

Type-II planets are brighter at smaller orbital distance 
so as to balance the larger stellar radiation, whereas these 
planets solidify on a shorter timescale because of the 
higher loss rate of water they experience. The sensitivity 
of the magma-ocean lifetime to the water amount and 
the orbital distance increases the possibility of placing 
constraints on the water endowment of terrestrial plan¬ 
ets based on the occurrence rate of molten planets for 
a given stellar age. Our model predicts that terrestrial 
planets remain molten over the main sequence of the host 
star in the type-II orbital region if the initial bulk content 
of water exceeds approximately 1 wt%. 

The composite spectrum is dominated by reflected 
light at the visible wavelengths and by thermal emis¬ 
sion at the near-IR wavelengths for type-II planets in the 
quasi-energy balance state. The atmospheric windows in 
the near-IR are fa vorable for d e tectin g thermal emission, 
as pointed out bv iLupu et al.l (|2014ll . The H-K color of 
the thermal emission is strikingly redder than that of the 
host star and the reflected light, whereas the solidifying 
type-I and type-II planets are typically as faint as planets 
having moderate surface temperature. Therefore, type- 
I and type-II planets could form a group distinct from 
the observed brown dwarfs and the planets in the solar 
system in the color-magnitude diagram. Evolution cal¬ 
culations including cloud formation will be required in 
order to make more detailed predictions regarding color 
and appearance. 

The K and L bands will be favorable for future di¬ 
rect imaging because the planet-to-star contrast of the 
type-II planets is approximately 10“^ — 10“®. In partic¬ 
ular, the search around Ocr, which is close to the inner 
edge of the habitable zone, will have a higher detection 
probability in terms of the duration of the molten state. 
For both the type-I and type-II planets, in visible atmo¬ 
spheric windows, the contrast of the thermal emission 
drops below 10“^® in less than 10® years, whereas that 
of the reflected stellar light remains 10“^®. Thus, the 
visible reflected light from molten planets will also be a 
promising target for direct imaging with ground-based 
high-contrast instruments and for second-Earths survey 
by space-based telescopes. 
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APPENDIX 

RADIATIVE TRANSEER 


We use two-stream approximation to solve radiative transfer equation in a plane-parallel atmosphere, with diffusivity 
factor rdf of 


dF'l 

rd{-^=F^ - J^{t) 

-r^,^=Fj-Mr). ( 1 ) 

Hereafter, the wavenumber subscript v is suppressed to simplify the notation. We assume that scattering is isotropic, 
and treat stellar and planetary radiation separately. Eor the planetary flux (eJj and , 

J(r) = Jpi(r) = (Uji + Fj,) + (1 - u;o)nB (T(t)) , (2) 

where ujq is the single scattering albedo, and B (T(t)) is the Planck function with temperature at t. Boundary 
conditions at the top and bottom of the atmosphere are given by 

Fp\, = 0 and Fp\i, = (l-Ag)7 rR(r.)+AgFp\j„ (3) 

where Ag is the ground albedo and Tg is the surface temperature. The subscripts t and b denote the top and the 
bottom of the atmosphere, respectively. 

Eor the stellar diffusive flux (fJ and Fj[) , 


J(r) = J.i(r) = (fI + Fi) + 


(4) 


where Dq is the stellar flux incident to the top atmospheric layer at optical depth of with zenith angle 6q (^q = cos 9q). 
The direct downward stellar flux is given by 


D4(r) = 1^01^406"/^“. (5) 

We assume that the stellar spectrum profile incident to the top of the atmosphere is given by that of a blackbody 
radiation with an effective stellar temperature T*. It gives the boundary condition as follows. 


F'*' — n 

Cglt ~ Cl, 


F>p, — Fr:-, 


aTX 


and fJ_|3 = Ag ( Fi k + D. 


sl,b 


sl.b 


( 6 ) 


where Fq is the total stellar energy flux that the planet receives at its orbital position. 

Wavenumber region is separated into two parts to calculate the radiative transfer equation efficiently and accurately. 
In the wavenumber region less than a certain wavenumber t'a, we assume that the scattering process by gaseous 



Fig. 11.— Spectrum of the thermal emission from the top of the steam atmosphere obtained using our model and that reported by 
IGoldblatt et all 1120131 ^ for a surface pressure of 123 bar and a surface temperature of 600 K. 
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Fi g. 12.— Planetary rad iation and Bond albedos for a steam atmosphere with a surface pressure of 260 bar. IKopparapu et ^ l|2Q13l') 
and IGoldblatt et al.l l|2013l^ obtained the Bond albedos by averaging the stellar flux over six zenith angles, whereas we calculated the Bond 
albedos using a stellar zenith angle of cos“^(l/\/3) (red closed circles and solid line). The open circles and dotted line represent the Bond 
albedos obtained by assuming that the stellar radiation is perfectly absorbed at wavenumbers beyond 30,000 cm~^. 


molecules can be neglected. The net upward planetary and stellar fluxes are calculated using the integrated form of 
Eq. O, using Eqs. (ED, O and (0, and the boundary conditions ([SD and (HD, 


FpVr) =7ri?(r(r)) + ^ ^ , 

F^r) =7rB(T(T)) - tt T g-r-ddt-Tl 

Jo dt 


-rdf(Tb-T) 


(7) 


and 


Esi(r)=0, 


( 8 ) 

given the fact that the diffusive downward stellar flux is equal to 0 in a non-scattering atmosphere. We evaluate the 
integrals in Eqs ED by considering the contributions from the region where the optical distance from the midpoint of 
each layer is smaller than 10. 

At sufficiently large wavenumber, both the absorption and scattering processes are taken into account. T he radiative 
flux in each layer is calculated using a formal solution of the radiative transfer equation m similarly to iToon et’ahl 
(|1989ll , but is solved iteratively by repeating upward and downward integrations to satisfy the boundary conditions at 
the top and the surface. 

In this paper, we calculate the radiative transfer equations for both the planetary and stellar radiation in a wavenum¬ 
ber range from 0 to 30,000 cm“^. The wavenumber z/^, at which the way to calculate the radiative transfer equation 
is switched from one to another, is set to 5,000 cm“^. Beyond 30,000 cm“^, only the scattering process is considered 
using the Rayleigh scattering cross-section at 30,000 cm“^. This leads to underestimate the attenuations by the ab¬ 
sorption and scattering processes, and overestimate the Bond albedo (see Eig. [TTD. We define a spectral albedo a for 
the incident stellar radiation as the ratio of the diffusive upward radiation to the direct downward radiation at the 
top of the atmosphere, i.e., Ej(Tt)/I?gj The Bond albedo for stellar radiation with an effective temperature of T* is 
obtained as a weighted mean of the spectral albedo as follows, 


Qn — 


J aB{T^)dv 


aTX 


(9) 


We compare the thermal emission spectr a of our model w i th tha t reported bv iGoldblatt et all (1201311 (Eig. [TTD. 
model spectra accord closely with that bviGoldblatt et aP (l2013f l. The thermal radiation calculated with our model 
is lower than that bv iGoldblatt et ahl (1201311 at the wavelengths between 3 and 5 ^m, due to the difference in the 
continuum absorption models of water vapor (see below). 

We also compare the planetary radiation and Bond albedos calculated with our radiative transfer model with those 
by previous studies to check the validity of our model (Eig. IT^. The radiation l imit o f our model is 280 W/m^, and 
is in close agreement with the value of 282 W/m^ obtained bv IGoldblatt et HI (l2013f l with line-by-line calculations, 
with a difference of less than 1%. As the surface temperature increases, the planetary radiation starts to rise for 
all the three models due to the t emper ature increase in the upper atmo s phere . In this phase, our results agree well 
with those by iKopparapu et al.l (1201311 , rather than by IGoldblatt et al.1 (1201311. proba.bly be cause of the differences 
in the treatment of the continuum absorption of water vapor. TPaynter fc Ramaswamvl (1201111 indica ted that the self 
continuum absorption of water vapor is underestimated in a MT_CKD 2.4 model, which was used bv iGoldblatt et ^ 
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(|2013fl . especially in a wavenumber range from 2000 to 3000 cm~ ^, compared with the MT_CKD 2.5 model used in 
our model and a BPS formalism used by iKopparapu et ahl (|2013|1. This like ly explains why the planetary radiation 
starts to increase at the lower surface temperature in iGoldblatt et al.l (|2013ll . The Bond albedos with our model are 
larger than those with the previous studies by about 0.05. This is probably because we underestimate the absorption 
and scattering processes at the wavenumbers larger than 30,000 cm“^. 
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